Interband superconductivity: contrasts between BCS and Eliashberg theory 
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The newly discovered iron pnictide superconductors apparently present an unusual case of interband-channel 
pairing superconductivity. Here we show that, in the limit where the pairing occurs within the interband channel, 
several surprising effects occur quite naturally and generally: different density-of-states on the two bands lead 
to several unusual properties, including a gap ratio which behaves inversely to the ratio of density-of-states; 
the weak-coupling limit of the Eliashberg and the BCS theory, commonly taken as equivalent, in fact predict 
qualitatively different dependence of the Ai/A 2 and A/T c ratios on coupling constants. We show analytically 
that these effects follow directly from the interband character of superconductivity. Our results show that in the 
interband-only pairing model the maximal gap ratio is y/Nj/Ni as strong-coupling effects act only to reduce 
this ratio. This suggests that if the large experimentally reported gap ratios (up to a factor 2) are correct, the 
pairing mechanism must include more intraband interaction than is usually assumed. 

PACS numbers: 74.20.Rp, 76.60.-k, 74.25.Nf, 71.55.-i 



Athough first proposed 50 years ago, multiband supercon- 
ductivity where the order parameter is different in different 
bands had not attracted much interest until 2001 when MgB2 
was found to be a two-band suprconductor. MgB2 represents a 
partucular case where one "leading" band enjoys the strongest 
pairing interactions, while the interband pairing interaction, 
as well as the intraband pairing in the other band, are weak. 
There is growing evidence that the newly disovered super- 
conducting ferropnictides represent another limiting case: the 
pairing interaction is predominantly interband, while the in- 
traband pairing in both bands is weak. This leads to a number 
of interesting and qualitatively new effects, including the fact 
that a repulsive interband interaction is nearly as effective in 
creating superconductivity as an attractive one. 

In this paper we will show another surprising feature of the 
two-band "interband" superconductivity (meaning supercon- 
ductivity induced predominantly by interband interactions): 
entirely counterintuitively, the BCS theory for such supercon- 
ductors is not the weak coupling limit of the Eliashberg theory, 
and the difference is not only quantitative but qualitative. This 
fact holds for either repulsive (as, presumably, in pnictides) or 
attractive interactions. 

Specifically, we will concentrate on the dependence of the 
superconducting gaps in the two bands on the ratio of the den- 
sities of states and the magnitude of the superconducting cou- 
pling. We will show that the gap ratio is always smaller in the 
Eliashberg theory than in the BCS theory, the deviation grows 
with coupling strength and with temperature, and is largest 
just below T c . 

Let us start with the BCS equations 1 . For a two band 
interband-only case, with gap parameters given on the two 



bands as Ai and A2, the BCS gap equations take the form 
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where E( j is the usual quasiparticle energy in band i given by 

£,•* — /j) 2 + A 2 , the normal state electron energy is e,-.*, /j 
is the chemical potential, and V is the interband interaction 
causing the superconductivity. V can be either attractive (> 
in this convention) or repulsive (as presumably in the pnic- 
tides), but for the rest of the paper the sign does not matter. 
For simplicity we will use V > and A > 0, keeping in mind 
that for pnictides all the results apply by substituting A by | A| . 
The BCS theory assumes V to be constant up to the cut-off en- 
ergy co 6 . Following the BCS prescription, we can convert the 
momentum sums to energy integrals up to a cut-off energy co c 
and assume Fermi-level density-of-states (DOS) Ni and N2. 
Near T c these equations can be linearized giving 



Ai = A 2 Xi 2 log(1.136G) c /r c ) 
A 2 = Ai^2ilog(1.136a> e /r e ), 



(2) 



where X\ 2 = N2V, the dimensionless coupling constant, with 
a similar expression for X 2 \ ■ These equatio ns readily yield 
h e ff = VX17X21 and a = A2/A1 = y/N\/N2- This result has 
been obtained before 2 ' 3 . Similarly, at T — in the weak- 
coupling limit 



Ai = A2A,i2sinh 1 ((o c /A2) 
A 2 = AiA.2isinh _1 (co c /Ai) 



(3) 
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Obviously, for X e /f — > we have T c — > and the relation 
A2/A1 = \/N{/N2 should hold. The same is not true for 

X eff > 0. 

First principle calculations suggest for the pnictides P = 
N2/N1 < 1.4, corresponding to the gap ratio a < 1.2. Experi- 
mental estimates for the gaps differ wildly, yielding gap ratios 
ranging from 1.3 to 3.4. Since the goal of this paper is to ad- 
dress the effect of the density of states difference on the gap 
ratio, we will use an intermediate number 4 a = 1 .6 (p = 2.6). 
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term can also be worked out and reads 
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FIG. 1 : (color online) The ratio of the gap functions in an interband 
pairing case, as a function of X e rr, for the BCS (dashed line) and 
Eliashberg Einstein (line) spectrum and spin fluctuation (triangle) 
spectrum cases. The dotted line represents numerical Eliashberg Ein- 
stein spectrum results in which the mass renormalization parameter 
has been artificially taken as 1, showing that the difference between 
BCS and Eliashberg is mainly a mass renormalization effect. Inset: 
analytic approximations to numerical results: diamonds are BCS Eq. 
[5] circles are Eliashberg Eq. \TE\ 



The fact that the band with the larger DOS ends up with 
a smaller gap is a somewhat counterintuitive result. This is 
a direct result of the interband-only pairing - the pairing am- 
plitude on one band is generated by the DOS on the other. 
Numerical solution of Eqs. [2] at T = (Fig. [T] gives, as ex- 
pected, a = = 1.61 at % e ff — > 0. As a function of X e /f 
it increases linearly, reaching s» 2.3 at X e ff « 1.6 (note that 
as shown below, it will ultimately saturate at P = 2.6 in the 
superstrong limit ). This increase can be easily explained. 

Let us define a = p 1 / 2 , so that x -C 1 at X -C 1 , and substi- 
tute sinh _1 (co e /A) — > log(2co c /A). A few lines of algebra then 



lead to 
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This result was also obtained by Bang 3 . The quadratic in X 
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^//(41ogp 



log 2 p) 



32 



(5) 



As Fig. Q] shows, this expression describes the numerical so- 
lution at small X very well. Although not apparent from the 
plots, the A1/A2 ratio will saturate at large X, as shown by 
Bang and Choi 3 and can also be seen from Eq. 3, since 

Ai = A 2 A.i2sinh _1 (co c /A 2 ) -> A,i 2 co 6 - forA 2 > co 6 (6) 

Similarly, in this limit A 2 = A, 2 i<a c so that Ai/A 2 = = 
N2/N1. All these BCS results, however are opposite to a 
known analytical result 5 that in the superstrong (Eliashberg) 
limit X 3> 1 the gap ratio a — > 1 independent of p. Let us 
now move to the strong-coupling limit, given by Eliashberg 6 
theory. 

In this theory, the BCS gap function Ay is replaced by 
a complex, energy-dependent quantity Ao(co), which must 
be determined along with a mass renormalization parameter 
Z(co). One commonly formulates the equations in terms of 
<j)(co) = Z(co)A(co), and these equations can be solved either 
on the real frequency axis or the imaginary axis (using Mat- 
subara frequencies). These equations are formulated in a two- 
band interband pairing case on the imaginary axis as follows 
(some of the notation is repeated from 7 ): 
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Here the kernel Ki 2 is given by 



Ku(i(Om - = 2 / 
Jo 



n 2 + (co„-co„) 2 



This B12 represents the electron-boson coupling function 
which supplants the pairing potential used in BCS theory, 
and there is an exactly analogous equation for band 2. Here 

fl U (fi)/2?2i(a)=tf 2 /tfi=P. 

First we assume a simple Einstein-type electron-boson cou- 
pling function. Numerical solution of the Eliashberg equa- 
tions dHJ finds that the ratio of the gaps decreases with X, 
opposite to the BCS prediction that the ratio of the gaps in- 
creases with increasing coupling. This can be understood an- 
alytically as well. 

First of all, we observe that neglecting the mass renormal- 
ization by setting Z = 1 in Eq|7]appears to be very close to the 
BCS solution (in fact, deviation from the lowest-order approx- 
imation of Eq. |4]is mainly due to the increasing difference be- 
tween sinh -1 (co c /A) and log(2co c /A)). Let us now work out 
the effect of the mass renormalization. 

Assum ing an Einstein spectrum with the frequency i2, at 
T=0 Eqs. |718lreduce to 



Ai(<b)Zi(<b) = 



/: 
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and 
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Zi(oo) = l + —X n Q. 2 [ da 

2(0 (fi 2 + (C0-Co') 2 )( X /co' 2 +A2(co)) 



with a similar equation for A2 and Z2. In the popular "square- 
well" approximation 8,9 the equations become 

A.mo = whcoD £^ x (9) 



A 2 (co ) 
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Zi(co) = 1 + — X n / <We(£2- Ico-co'l) x 
2co _/_<*, 
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(10) 



which may be readily integrated to yield the following renor- 
malization behavior for Z(co): 

Zi(co) = 1 +A,i 2 forco < i2 (11) 
= l+A,i 2 n/cofori2<co<2n (12) 
= l+Xi 2 /2forco>2n (13) 

This mass renormalization behavior can then be incorporated 
in the previous BCS equations yielding a natural result: 

Ai(l +X12) = A 2 >.i2sinh- 1 (co c /A 2 ) (14) 
A 2 (l+M = Ai^sinh-^cojAi), (15) 

reducing to Eq. |3]with A,i 2 — > Xi 2 /(1 + A,i 2 ),A, 2 i — > A, 2 i/(1 + 
A, 2 i). Thus, in the linear order in X. 



X21 — X 



12 



(16) 



The last term is negative and always larger than the previous 
one (independent of p). Thus, the net effect is always oppo- 
site to what the BCS theory predicts. We have plotted up the 
above analytic approximation in Figure 1 (solid line in inset) 
and find good agreement for X e /f < 0.4, showing that the mass 
renormalization is responsible for the lessening of the gap ra- 
tios with increasing coupling in Eliashberg theory. This result 
might in hindsight have been expected given that the Fermi 
surface with the larger gap at weak-coupling can be expected 
to have larger self-energy interactions in Eliashberg theory, re- 
ducing the gap anisotropy. This result is also consistent with 
the superstrong coupling limit of equal gaps, as mentioned 
previously. 

Interestingly, this strong coupling effect remains operative 
at all temperatures up to T c , while the previous term in Eq. 
[I6]vanishes at T c . Therefore (cf. Fig|2]i the actual gap ratio is 
even closer to 1 near T c than at T = 0. 

Finally, we note that the above Eliashberg results were ob- 
tained using an Einstein spectral function for simplicity, but 
as indicated on the plot the use of a typical spin-fluctuation 
spectrum [~ coi2/(<B 2 +£2 2 )] does not alter the results. 



Another interesting observation to be made concerns the 
A(0)/T c ratios predicted by BCS and Eliashberg theory. In 
the conventional weak-coupling one-band BCS theory this ra- 
tio does not depend on X. This is no longer the case in the 
two-band BCS with the interband coupling only. In the low- 
est order the reduced gaps are simply Ai(0)/T ( . = 1.760 1 / 4 , 
A 2 (0)/r c . = 1.76p-'/ 4 . The next order can be worked out us- 
ing Eq.|5] 

m = L76p ./. (1+ xlM_i2££ ) (1V) 

4> - 1.76 P ^(l-* 4l ° gP + '° 82P ) ,18, 

This is confirmed by numerical calculations (Fig. [2}: the 
smaller gap ratio decreases with X, while the other gap in- 
creases. Since the Eliashberg equation makes the gaps closer 




FIG. 2: (color online) A(T = 0) /T c ratios are shown, as a function of 
the overall coupling constant, for the BCS (solid black line and gray 
dashed line) and Eliashberg (dot-dashed line and double-dot-dashed 
line) cases. 

with increased coupling, this odd behavior does not show up: 
both reduced gaps grow with A.. 

For completeness, we also show in Fig. 3 the behavior 
(in Eliashberg theory) of the reduced gaps as a function of 
the DOS ratio N%/N\ =X\2/X2\- As might be expected, as 
the DOS ratio becomes very small the gap ratios move apart 
appreciably. Interestingly, T c (shown in the right panel) is not 
constant as it would be in a weak-coupling regime, but varies 
significantly for coupling constant ratios far from 1 . This is 
a result of the use of comparatively large coupling constants 
on one band when the other coupling constant is small, so that 
T c suppression due to thermal excitation of real phonons (an 
effect not present in the BCS formalism) is stronger. 

To conclude, in this work we have shown for the interband- 
only pairing the two-band superconductivity is qualitatively 
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incorrectly described by the BCS formalism even for the weak 
coupling limit. BCS and Eliashberg theory predict qualita- 
tively different behavior (as a function of coupling constant) 
for such basic characteristics as the gap ratio a = A1/A2, as 
well as for the reduced gaps A/T c . In particular, the sign 
of da/dX changes from BCS to Eliashberg theory. We 
have found this result analytically and numerically, by solv- 
ing Eliashberg equations for model spectra. This finding is 
relevant to the superconducting pnictides where the interband- 
pairing regime is believed to be realized. 



FIG. 3: (color online), (left) The behavior of the Eliashberg A(0)/r c 
ratios as a function of the ratio of coupling constants, (right) The 
behavior of T e in this case. For both cases A, e ff is fixed at 1. 
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